Mixtures of Bose gases under rotation 
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We examine the rotational properties of a mixture of two Bose gases. Considering the limit of 
weak interactions between the atoms, we investigate the behavior of the system under a fixed angular 
momentum. We demonstrate a number of exact results in this many-body system. 
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One of the many interesting aspects of the field of cold 
atoms is that one may create mixtures of different species. 
The equilibrium density distribution of the atoms is an 
interesting problem by itself, since the different compo- 
nents may coexist, or separate, depending on the value 
of the coupling constants between the atoms of the same 
and of the different species. If this system rotates, the 
problem becomes even more interesting. In this case, the 
state of lowest energy may involve rotation of either one 
of the components, or rotation of all the components. 
Actually, the first vortex state in cold gases of atoms was 
observed experimentally in a two-component system 
following the theoretical suggestion of Ref . [2] . More re- 
cently, vortices have also been created and observed in 
spinor Bose- Einstein condensates [3, 0|. Theoretically, 
there have been several studies of this problem [!, 0, 0] , 
mostly in the case where the number of vortices is rel- 
atively large. Kasamatsu, Tsubota, and Ueda have also 
given a review of the work that has been done on this 
problem [§]. 

In this Letter, we consider a rotating two-component 
Bose gas in the limit of weak interactions and slow ro- 
tation, where the number of vortices is of order unity. 
Surprisingly, a number of exact analytical results exist 
for the energy of this system. The corresponding many- 
body wavefunction also has a relatively simple structure. 

We assume equal masses M for the two components, 
and a harmonic trapping potential V t — M(u> p + 
u> 2 z )/2, with p 2 — x 2 + y 2 . The trapping frequency 
lo z along the axis of rotation is assumed to be much 
higher than u. In addition, we consider weak atom- 
atom interactions, much smaller than the oscillator en- 
ergy Hlo, and work within the subspace of states of the 
lowest Landau level. The motion of the atoms is thus 
frozen along the axis of rotation and our problem be- 
comes quasi-two-dimensional The relevant eigen- 
states are $ m (p, 0)<po(z), where <fr m (p,8) are the lowest- 
Landau-level eigenfunctions of the two-dimensional os- 
cillator with angular momentum mh, and (fo(z) is the 
lowest harmonic oscillator eigenstate along the z axis. 

The assumption of weak interactions also excludes the 
possibility of phase separation in the absence of rotation 
(lfj | , since the atoms of both species reside in the lowest 
state $o.o( r ) = &o(p, 8)<Po(z), while the depletion of the 
condensate due to the interaction may be treated pertur- 
batively. 

We label the two (distinguishable) components of 



the gas as A and B. In what follows the atom- 
atom interaction is assumed to be a contact potential 
of equal scattering lengths for collisions between the 
same species and the different ones, ciaa = &bb = 
clab = o< The interaction energy is measured in units 
of v a = Uq J |$ ,o(r)| 4 d 3 r = (2/iT) 1 / 2 Tu J ja/a z , where 
U = 4:nh 2 a/M, and a Q = (S/Afa) 1 ' 2 , a z = {h/Muj^ 1 ' 2 
are the oscillator lengths on the plane of rotation and 
perpendicular to it. 

If Na and Nb denote the number of atoms in each 
component, we examine the behavior of this system for 
a fixed amount of L units of angular momentum, with 
< L < N max , where A max = max(N A , Nb)- We use 
both numerical diagonalization of the many-body Hamil- 
tonian for small systems, as well as the mean-field ap- 
proximation. Remarkably, as we explain in detail below, 
there is a number of exact results in this range of angular 
momenta. 

More specifically, when < L < N m i n , where N m i n = 
min(A r J 4, Nb), using exact diagonalization of the many- 
body Hamiltonian, we find that the interaction energy of 
the lowest-energy state has a parabolic dependence on L 
in this range, 



£o(L)/v = ^N(N 



1 



Inl + \l{l 



1 



(1) 



with N = Na + Nb- In addition, the lowest-energy state 
consists only of the single-particle states of the harmonic 
oscillator with m = and m = 1 . The occupancy of the 
m = 1 state of each component is given by 

N B - L + 1 
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while (N A ) m =o 
{N B ) m =i- 
As L becomes larger than A m i n , 
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there is a phase 



transition. For A min < L < N max , the single-particle 
states that constitute the many-body state are no longer 
only the ones with m — and m = 1 (as in the case 
< L < N min ), and in addition the interaction energy 
varies linearly with L, 
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FIG. 1: The interaction energy that results from numerical 
diagonalization of the Hamiltonian of a mixture of two Bose 
gases, with Na = 4 and Nb = 12 (lower curve, marked by 
"+"), as well as Na = 16, and Nb = (higher curve, marked 
by "x"), as function of the angular momentum L, for < L < 
12. 
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FIG. 2: The occupancy of the single-particle states with 
m — 0,1 and 2, as function of the angular momentum 
< L < 12, that results from numerical diagonalization of 
the Hamiltonian of a mixture of two Bose gases, with Na = 4 
and Nb ~ 12. The upper panel refers to species A, and the 
lower one to species B. 



The lower curve in Fig. 1 shows the interaction energy 
of a system with Na = 4 and N B = 12, for < L < 
12. For < L < 4 the energy is parabolic, and for 
4 < L < 12, it is linear. These are exact results, within 
numerical accuracy. The higher curve is the interaction 
energy of a single-component system of N — 16 atoms. 
It is known that in this case, the interaction energy is 
exactly linear for 2 < L < N = 16 This line is 

parallel to the line which gives the interaction energy 
of the system with Na — 4 and Nb — 12 for Na — 
4 < L < Nb = 12. Figure 2 shows the occupancy of 
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FIG. 3: The conditional probability distribution of a mixture 
of two Bose gases, with Na = 4 (higher row), and Nb = 12 
(lower row). Each graph extends between — 2.4ao and 2.4ao. 
The reference point is located at (x,y) = (ao,0) in the higher 
graph. The angular momentum L increases from left to right, 
L — 2,3, 4(= Na), 5, 6, 8, 10, 11, and 12(= N B ). 



the single-particle states, that result from the numerical 
diagonalization of the Hamiltonian. 

The physical picture that emerges from these calcula- 
tions is intriguing: as L increases, a vortex state enters 
the component with the smaller population from infinity 
and ends up at the center of the trap when L — N m [ n . In 
addition, another vortex state enters the component with 
the larger population from the opposite side of the trap, 
reaching a minimum distance from the center of the trap 
for L « N A /2 (this estimate is valid if 1 < N A < N B ), 
and then returns to infinity when L = A m j n . This mini- 
mum distance is « 2(Nb/Na)cio. For N min < L < iV max , 
the vortex in the cloud with the smaller population (that 
is located at the center of the trap when L = N min ) moves 
outwards, ending up at infinity when L — iV max , while a 
vortex in the other component moves inwards again, end- 
ing up at the center of the trap when L = N^x. Figure 
3 shows clearly these effects via the conditional proba- 
bility distributions, for Na = 4, and Nb = 12. We note 
that in the range < L < Na = 4, the (distant) vortex 
in the large component (lower row, species B) is too far 
away from the center of the cloud to be visible, because 
of the exponential drop of the density. The plots in Fig. 3 
(and Fig. 5) are not very sensitive to the total number of 
atoms N = Na + Nb , and resemble the behavior of the 
system in the thermodynamic limit of large N. 

In the case of equal populations, Na — Nb, the 
parabolic expression for the interaction energy, Eq. |T]). 
holds all the way between < L < Na — Nb- Figure 4 
shows the interaction energy and Fig. 5 the occupancies 
of the single-particle states, which vary linearly with L. 
The corresponding physical picture is quite different in 
this case, as shown in Fig. 6. The system is now sym- 
metric with respect to the two components, and a vortex 
state enters each of the components (from opposite sides) . 
These vortices reach a minimum distance from the center 
of the trap equal to do, when L — Na = Nb- 

The simplicity of the system that we have studied al- 
lows one to get some relatively simple analytical results, 
which we present below. As we saw earlier, when L = Na 
or L = Nb, there is a unit vortex state in species A or 
B, while the other species is in the lowest oscillator state 
with m = 0. Since (at least to leading order and next to 
leading order) only the states with m = and m = 1 are 
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FIG. 4: The interaction energy that results from numerical 
diagonalization of the Hamiltonian of a mixture of two Bose 
gases, with Na = Nb = 8, as function of the angular momen- 
tum L, for < L < 8. 
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FIG. 5: The occupancy of the single-particle states with m = 
and m = 1, as function of the angular momentum < 
L < 8, that results from numerical diagonalization of the 
Hamiltonian of a mixture of two Bose gases, with Na = Nb = 
8. 

occupied, the Fock states are of the general form (if, for 
example, L = N B ) 

|n) = |0^-™, 1") (g) |0", l^-"). (5) 

Expressing the eigenstates of the interaction V as |$) = 
5^ n (— l)"/ n |n), the eigenvalue equation takes the form 

Vn,nfn ~ Vn,n-lfn-l ~ Vn,n+lfn+l = £ fn , (6) 

where V n , m are the matrix elements of the interaction 
between the above states. Remarkably, if Na = Nb = 
N/2, then 

V n , n - V n>n - X - V n , n+ i = 5N(N - 2)u /16, (7) 

which implies that in this case, the lowest eigenen- 
ergy is £q — 5N(N — 2)uo/16, in agreement with 
Eq. |T]). The corresponding eigenfunction is simply 

N) =£„(-!)»• 

In the case Na ^ Nb, with, e.g., L = Nb, n is of order 
unity, and therefore the interaction may be written as 

V/v = \n a {N a - 1) + \n b {N b - 1) + \n a N b + 

+-(N A + N B )ala + -\/N A N B (a b\ + aJSi),(8) 



FIG. 6: The conditional probability distribution of a mixture 
of two Bose gases, with N A = Nb = 8. The two rows refer to 
the two different species. Each graph extends between — 2.4ao 
and 2.4a . The reference point is located at (x,y) = (a ,0) 
in the lower graph. The angular momentum L increases from 
left to right, L = 2,3, ... ,8. 

where a m (aj n ) and 6 m (6j„) are annihilation (creation) op- 
erators of the species A and B with angular momentum 
rrih. The above expression for V can be diagonalized 
with a Bogoliubov transformation, 

V/v = \\N A - N b \(2&& + 1) + ^N A {2N A - 3) + 

+\n b {Nb-2) + ^N a N b , (9) 

where a^a is a number operator. When N A = N B = 
N/2, the lowest eigenenergy is 5N(N — 2)/16, in agree- 
ment with Eq. ([I]) . When N A ^ N B , the lowest eigenen- 
ergy is 

£o/vq = \\N A - N B \ + \n a {2N a - 3) + 

+\n b {N b -2) + ^N a N b . (10) 

The above expression agrees exactly with Eq. ([I]) when 
N A < Nb = L, and to leading order in TV when N A > 
N B =L. 

In addition, according to Eq. the excitation ener- 
gies are equally spaced, separated by — N B \v /2 + 
0(vq). Therefore, one very important difference between 
the case N A = Nb and Na ^ Nb is that in the first 
case there are low-lying excited states, with an energy 
separation of order vq , while in the second [where in gen- 
eral N A — Nb ~ O(N)], the low-lying excited states are 
separated from the lowest state by an energy of order 
Nv . 

Let us now turn to the mean-field description of this 
system, for < L < N min . We consider the following 
order parameters for the two species (restricting ourselves 
to the states with m — and m = 1 only), 

* A = (co$o + ci$i)^ (z), * B = W)*o + di$x)<Po(z), 

(11) 

where cq, c\, do, d\ are variational parameters. Given 
the order parameters, the many-body state is $mf = 
Ililfi &A{ r i) Ylf=i ^ > s( r i)- The normalization for each 
species implies that |co| 2 + | ci | 2 = 1, and \do\ 2 + \di\ 2 = 
1, while the condition for the angular momentum gives 
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I ci | 2 + Nb\<1i\ 2 = L. The interaction energy is 
1 



£mf = -N A (N A -1)U J \<S> A \U 3 r + 
+ \n b (N b -1)U q J \<f B \ 4 d 3 r + 
+ N A N B UoJ \3> A \ 2 \$ B \ 2 d 3 r, (12) 
or, to leading order in N, 
£mf 1 / N A 



N 2 v 2 V N 



= 7rKf icor + dcir + sMicin + 



. N A N B (. |2 , , ,2 . 1 |2|j |2 , 1| |2|j |2 , 

+ — (|co| |do| +2l c °l ' dl + 2' Cl ' ' 

+ ^|ci| 2 Mi| 2 -|co|| C i||dolMi|),(13) 

where the phases of the variational coefficients have been 
chosen so as to minimize £mf- 

When N A = N B , then |c | 2 = |d | 2 , and also |c x | 2 = 
|di| 2 , which implies that |ci| 2 = |di| 2 = L/N = I. There- 
fore, one finds that £mf/N 2 = (2 — 2l + l 2 )v /4:, in agree- 
ment (to leading order in N) with the result of exact 
diagonalization, Eq. (p}. When Na ^ Nb, minimization 
of the energy with respect to one of the four (free) varia- 
tional parameters (the other three are then fixed by the 
three constraints) gives a result that agrees to leading 
order in N with that of numerical diagonalization. 

The fact that for < L < AT m j n the lowest many-body 
state consists of only the m — and m — 1 single-particle 
states is remarkable. To get some insight into this result, 
we consider the two coupled Gross-Pitaevskii equations, 
which describe the order parameters ^>a and ^b- If MA 
and [Ib is the chemical potential of each component, then 



ft 2 v 2 

' 2M 



V t + Uo\y B \ 



U \^ A \ = MA*A, 



h 2 v 2 

2M 



Vt + Uo\^ A \ 2 ) ^b + U \^b\ 2 ^b = 

(14) 



For a large population imbalance, where, for example, 
1 < N A < Nb, for < L < N min = N A , most of 
the angular momentum is carried by the species A (with 
the smaller population). As mentioned earlier, in the 
range < L < N m \ n , although there is also a vortex 
state in species B, this is far away from the center of 
the cloud. As a result, the order parameter of species B 
is essentially the Gaussian state, with the corresponding 
density being \$> B {p, z)\ 2 «ns(0, 0)e~ p / a o- z / a * ; where 
ns (0,0) is the density of species B at the center of the 
trap, i.e., at p — and 2 = 0. 

This component acts as an external potential on 
species A. Thus, the total "effective" potential acting 
on species A is (expanding the function e~ p / a °), 



M 



\ ;,,(,..( 1 1 * + cw(o,o) ^i - ^ + ^ ) .( 



for distances close to the center of the cloud. One may 
argue that in this self-consistent analysis, the quadratic 
term in the expansion changes the effective trap fre- 
quency, while the quartic term acts as an anharmonic 
potential. We argue that this anharmonic term is re- 
sponsible for the fact that only the states with m = and 
m = 1 are occupied [13]. Actually, this is more or less 
how Dalibard et al. investigated the problem of multiple 
quantization of vortex states [l3|. In that case, it was an 
external laser beam that created an external, repulsive 
Gaussian potential, as opposed to the present problem, 
where this potential results from the interaction between 
the different components. 
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